Out-of-Equilibrium Kondo Effect: Response to Pulsed Fields 
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The current in response to a rectangular pulsed bias potential is calculated exactly for a special 
point in the parameter space of the nonequilibrium Kondo model. Our simple analytical solution 
shows the all essential features predicted by the non-crossing approximation, including a hierarchy of 
time scales for the rise, saturation, and fall-off of the current; current oscillations with a frequency of 
eV/h; and the instantaneous reversal of the fall-off current for certain pulse lengths. Rich interference 
patterns are found for a nonzero magnetic field (either dc or pulsed), with several underlying time 
scales. These features should be observable in ultra small quantum dots. 



PACS numbers: 72.15.Qm, 72.10.Fk, 73.50.Mx 

The recent observation of the Kondo effect in ultra 
small quantum dotsErlj has focused considerable atten- 
tion on the Kondo effect in mesoscopic systems. Acting 
as tunable Anderson impurities, quantum-dot devices of- 
fer an outstanding opportunity to test some, of the ba- 
sic notions of the Anderson impurity model.B For exam- 
ple, the crossover from the local-moment to the mixed- 
valent and empty-impurity regimes, El the tt/2 phase shift 
associated with the KondjCu,effectj3 r arid_the evolution of 
the Kondo effect with dcLTu and adi9iI3 bias potentials. 
Most recently, an intriguing scenario for the response of a 
quantum dot to a rectangula&,pulsed bias potential was 
put forward by Plihal et al£3 Using a time-dependent 
version of the non-crossing approximation, these authors 
predicted a hierarchy of time scales for the rise, satura- 
tion, and fall-off of the current through the dot, featuring 
a rich interplay between the Kondo temperature, the ap- 
plied voltage bias, and the temperature. For larger values 
of the pulsed bias eV, oscillations were predicted in the 
time-dependent current, with a characteristic frequency 
of lu = eV/h. These current oscillations were attributed 
to the sharp excitation frequency between the ±wo split 
Kondo peaks induced by the eV bias potential.EJ 

In this paper, we exploit an exactly solvable nonequi- 
librium Kondo modefiEJ to obtain a simple analytic de- 
scription of the Kondo-assisted current in response to 
rectangular pulsed fields. Our solution captures all the 
essential features reported by Plihal et aZ.J13 including 
the hierarchy of time scales for the rise, saturation, and 
fall-off of the current; the eV/h oscillations in the time- 
dependent current; and the instantaneous reversal of 
the current - for certain pulse lengths - after the pulse 
has ended. Furthermore, we are able to include also a 
nonzero magnetic held which itself can be pulsed, reveal- 
ing complex interference patterns between the bias and 
magnetic field. Most importantly, since our solution is 
analytic and exact, it provides a benchmark result for 
the response of a quantum dot to pulsed fields. 

We begin with a few general remarks on the solvable 



model used. An extension of the ToulouseEj and Emery- 
Kivelsonlij solutions, of the equilibrium Kondo problem 
to nonequilibrium,E3 this model is expected to correctly 
describe the strong-coupling regime of the nonequilib- 
rium Kondo, effect. -Ludeed, previous applications of the 
model to dda and adB bias potentials have shown all the 
qualitative features of Kondo-assisted tunneling: A zero- 
bias anomaly that splits in an applied magnetic field; 
Fermi-liquid characteristics in the low-T and low-V dif- 
ferential conductance; Side peaks in the differential con- 
ductance at eV = ±nhuj, for an ac drive of frequency 
to. At the same time, since the solvable point involves 
strong couplings, this model is incapable of describing 
weak-coupling features such as the logarithmic tempera- 
ture dependence of the conductance at elevated T. 

Introducing the actual physical system, it consists of 
two noninteracting leads of spin- 1/2 electrons, interact- 
ing via a spin-exchange coupling with a spin- 1/2 impurity 
moment f placed in between the two leads. The impurity 
moment can represent either an actual magnetic impu- 
rity or a quantum dot with a single unpaired electron. 
Each lead is subject to a separate time-dependent po- 
tential, such that a time-dependent voltage bias forms 
across the junction. Since the interaction with the impu- 
rity is local in space (s-wave scattering), one can reduce 
the conduction-electron degrees of freedom that couple to 
the impurity to one-dimensional fields ij} a(J {x), where a = 
R,L labels the lead (right or left) and a =f,J. specifies 
the spin orientations. In terms of these one-dimensional 
fields, coupling to the impurity takes place via the local 
spin densities s a/3 = \ YLa,o' ^L-(^^°-''^jW( )- 

As discussed in detail in Refs. |j and [hy, the exactly 
solvable model corresponds to a particular choice of the 
spin-exchange couplings, described by the Hamiltonian 

H = H - eV a {t)N a - n B g t H{t)T z 

a=L,R 
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Here the inter-lead longitudinal couplings have been set 
equal to zero, while the intra-lead longitudinal couplings 
take the value J z = 2-kTivf- The three transverse cou- 
plings Jf L , J^ R , and jf R = jj^ are arbitrary. N a and 
V a (t) in the second term of Eq. (|j) are the number opera- 
tor and time-dependent potential on lead a, respectively, 
while H(t) is a local magnetic field acting on the impurity 
spin. 

The basic approach we take is identical to that of 
Rcf. [H], where the same system was solved for ac drives. 
The first step is to convert the interacting nonequilibrium 
problem to a noninteracting one. To this end, we first de- 
compose the Hamiltonian into an unperturbed part and 
a perturbation, where the unperturbed part consists of 
Ho and the time-dependent voltage V a (t). The initial 
density matrix — characterizing the distribution of the 
unperturbed system before tunneling has been switched 
on — is equal to po — e~^"° /Trace {e _,8Wo }. 

The Hamiltonian, its unperturbed part, and the ini- 
tial density matrix are then convertfteLto quadratic forms 
via the following sequence of steps P'E3 (i) Bosonizing the 
fcrmion fields; (ii) Converting to new boson fields cor- 
responding to collective charge (c), spin (s), flavor (/), 
and spin-flavor (s/) excitations; (iii) Performing a canon- 
ical transformation; (iv) Refermionizing the new boson 
fields by introducing four new fermion fields, VvW with 
ji = c, s, /, sf; and (v) Representing the impurity spin 
(which has been mixed by the canonical transformation 
with the conduction-electron spin degrees of freedom) 
in terms of two Majorana fermions: a — —\f2T y and 
b = — \[2t x . For J z — 2ttUvf, one thus arrives at the 
following representation of the Hamiltonian: 

H' = E e *4 A* " * V c {t)N c + eV(t)N f 

l*=C,B,f,sf k 

+ifi B g i Hba + i — ^= V7^L k + i> a f,k)b 
2y/TTaL 
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Here J ± are equal to (J^ L ± J RR )/2; the energies et 
are equal to hvpk; N c and Nf are the number operators 
for the charge and flavor fermions, respectively; V c (t) and 
V(t) are equal to V L (t) + V R {t) and V R (t)-V L [t), respec- 
tively; a is an ultraviolet momentum cutoff, correspond- 
ing to a lattice spacing; and L is the effective size of the 
system (i.e., k is discretized according to k = 2ttti/L). 
Within this mapping, the Hamiltonian term Tia is re- 
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FIG. 1. The time-dependent current I(t) in response to 
an eV bias potential that is switched on at time t = 0, 
and then switched off again at time t = t. Here r = 5tx 
with tK = h/Ta, ksT is fixed at 0.1F a , and Idc(V) is the 
steady-state current for a dc voltage bias of eV [see Eq. (m)]. 
As eV is increased, oscillations develop in I(t), leading to a 
significant reduction in the rise time. The saturation time, 
however, is unaffected by the bias at the solvable point. In- 
set: Temperature dependence of I(t) for fixed eV/T a = 4 and 
i"/tic = 3. Full, long-dashed, dashed, and dotted lines cor- 
respond to knT/Ta = 0.1,0.5,1 and 2, respectively. With 
increasing T, the saturation time is reduced and the current 
oscillations are suppressed. 



and so po is transformed to e~ /3w '=»«/Trace{e _ ' 3w '= i «}. 

Having converted the problem to a noninteracting 
form, one can sum exactly all orders in the perturbation 
theory to obtain the.|Current. The solution features two 
basic energy scales,^ T a = [(Jf 77 ) 2 + ( J-) 2 ]/ AiraTiVF 
and Tb = ( J + ) 2 / AitafiVF, which determine the width of 
the different Majorana spectralrfunctions, and thus play 
the role of Kondo temperatures £3 The conventional one- 
channel Kondo effect is best described by the case where 
r a = Tb, for which only a_aingle Kondo scale emerges B 

Following Plihal et al.jEl we begin with the case of a 
pulsed bias potential and a zero magnetic field. Explic- 
itly, we assume that a bias potential of eV^ is switched 
on at time t — 0, and then switched off again at time 
t = t. Denoting (J± R ) 2 /AnahvF by Ti, the resulting 
time-dependent current I(t) flowing from right to left is 
conveniently expressed via the auxiliary function 



-[(2n+l)Trk B T+r a ]t/h 
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duced to the free kinetic-energy part of Eq. 
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For t < 0, I(t) is obviously zero. For < t < r it equals 
h{t) = I dc (V) - Im {e- ieVt ' h J{t; 0, V)} , (5) 
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while for r < t it is given by 

h(t) = lm{J(t-r; 0, V)} - Im [e~ ie VT ' h J{t; 0, V)} . 

(6) 

Here IdcO{) is the steady-state current for an eV dc volt- 
age bias:El 



dcV ; irh I V 2 2nk B T 
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The relevant time scales for the evolution of the current 
are easily read off from Eqs. (||) and (||). As proposed by 
Plihal et al.,c3 for < t < r the current has two compo- 
nents: a constant term equal to the new steady-state cur- 
rent, and a transient term that oscillates with frequency 
lj — eV/Ti. The latter component involves two natural 
times scales: the period of oscillations t v = 2nfi/eV , and 
the basic decay time td = h/(Y a + TtksT), which enters 
via J(t,0, V). For td <C t v , the current saturates before 
any oscillations develop. Hence the rise and saturation 
times are both of the order of td- On the other hand, for 
t v < td the current oscillates before saturating at Idc{V), 
resulting in a much shorter rise time. Specifically, the 
rise time is identified with the first instance at which the 
transient term vanishes. Writing J(t; 0, V) of Eq. (||) as 
| J(t; 0, V)|e^W, it is simple to show that < <p{t) < f . 
As the transient current vanishes each time eVt/h— ip(t) 
equals an integer multiple of 7r, this first happens at some 
time shorter than t v /A. Thus, while the saturation time 
still follows td, the rise time is bounded by t v /4 <C td- 
This transition in rise time from approximately td to less 
than t v /4 <C td is clearly seen in Fig. |l|, where the time- 
dependent current is plotted for increasing values of eV. 

Contrary to the current oscillations that may develop 
for < t < t, the fall off for t > t is described in Eq. (Q) 
by the difference of two non-oscillatory terms, each of 
which decays to zero. While the tail of the fall off is 
always characterized by the fall-off time td, the initial 
stages of the drop do not typically follow a single time 
scale. This is best seen for long pulse durations, r 3> td, 
when only the first term survives in Eq. (J6j) . To this 
end, consider Im{ J(t — t; 0, V)} with t = t and t v < td- 
Using the series representation of Eq. (Q) , an increasing 
number of terms contribute to lm{J'(0;0,V)} as eV is 
increased, typically up to n ~ eV/2itkBT. For t > t, 
each of these terms decays with a separate relaxation 
time, t d n) = h/[(2n + l)nk B T + T a ], such that the total 
signal contains a range of time scales from ~ t v up to td- 
While all time scales participate in the early stages of the 
fall off, the tail is governed by the longest time scale, td- 

For shorter pulse durations, r ~ t v < td, the fall-off 
current is actually quite sensitive to the precise value of 
r. This is demonstrated in Fig. ||, where several pulse 
durations are plotted for fixed eV/T a — 10. Although 
eV/T a = 10 may appear too large a bias for the applica- 
tion of the solvable point, the qualitative agreement with 
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FIG. 2. The current I(t) in response to a pulsed bias po- 
tential of magnitude eV — 10r a , for different pulse dura- 
tions r. Here k B T = 0.ir a and t v = 2?r?i/eV. I dc (V) is the 
steady-state current for an eV dc voltage bias. The asymp- 
totic fall-off current oscillates with sin(2nT /t v ) , resulting in 
the reversal of the current for r/t v = 3/4. 



the numerical calculations of Plihal et al. (see Rcf. [TJ, 
Fig. 3) is strikingly good. In fact, we are even able to 
reproduce the reversal of the fall-off current for certain 
pulse lengths. To understand the origin of this somewhat 
surprising effect, we go back to Eq. (o) with a sufficiently 
large t — r, such that each J function can be roughly 
approximated by the n = term in Eq. (0) (we assume 
T > 0). For a large voltage bias such that eV ^S> T a and 
eV ^> irksT (as is the case in Fig. ||), the fall-off current 
of Eq. (0) is asymptotically given by 
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Thus, since t d 3> t v and r ~ t v under the assumptions 
above, the asymptotic fall-off current is dominated by 
the last term in Eq. (JsJ> . Consequently, the sign of hit) 
oscillates with sin(27rr/t„), precisely as seen in Fig. 0. 

Although our solution clearly captures the essential 
findings of the non-crossing approximation, there are 
two major omissions to be noted. First, working with 
a Kondo impurity rather than an Anderson impurity, 
our model lacks the short charge-fluctuation time scale 
(to = h/Tdot hi the notations of Ref. |l3|), which governs 
the very early response of a quantum dot to the abrupt 
change in the applied voltage bias. Second, for < t < r, 
our solution does not show the expected decrease in satu- 
ration time for large voltage biasJiia This decrease, which 
stems frop, the dissipative lifetime induced by the bias 
potential,!!-] is not captured by the solvable point. On the 
other hand, our solution gives a very transparent picture 
for the role of a temperature. As the temperature is in- 
creased, td crosses over from Ti/Y a to h/irksT, which re- 
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FIG. 3. The time-dependent current for (a) a dc magnetic 
field H and a rectangular pulsed bias potential of magnitude 
eV — 4F a ; and (b) a dc bias potential eV = 2F a and a 
rectangular pulsed magnetic field of magnitude H. Upper 
panel: Solid, dashed, dotted, and long-dashed lines corre- 
spond to fiBgiH /T a = 0, 1, 2, and 4, respectively. Lower 
panel: Solid, dotted, and long-dashed lines correspond to 
fiBgiH/r a = 1,2, and 4, respectively. In both cases, the 
pulse duration is r = 5tk (t^ = fi/Y a ), and the temperature 
is equal to kBT/Y a = 0.1. Idc(V,H) marks the steady-state 
current for a dc voltage bias of eV and a dc magnetic field H. 



mains the only relevant time scale for ksT 3> eV. Thus, 
all response times for the current are determined by the 
temperature in this limit. 

An important advantage of the solvable point with 
r a = Tb is the ability to incorporate an arbitrary time- 
dependent magnetic field. Specifically, we find that the 
current I(t) for a general voltage bias eV(t) and an arbi- 
trary field H(t) is equal to \ [J+(f) + /-(f)], where I±(t) 
is the current for a zero magnetic field and an effective 
bias potential eV±(t) = eV(t) ± fj,ggjH(t). The latter 
bias has a simple physical interpretation. When an elec- 
tron tunnels by flipping the impurity spin, the effective 
potential barrier it sees is equal to eV±(t), depending on 
the initial impurity state, and the direction of tunneling. 

Using the above decomposition of I(t) into I±(t), we 
have computed the current for the two opposite cases of 
(i) a dc magnetic field H and a rectangular pulsed bias 
potential of magnitude eV, and (ii) a dc bias potential of 
eV and a rectangular pulsed magnetic field of magnitude 
H . In both cases, the pulse duration was taken to be r, 
with the corresponding field being equal to zero for t < 
and t > t. The calculation of I±(t) in each of these cases 
requires the generalization of Eqs. (||)-((6|) to a pulsed 
bias potential of the form V(t) = V2 for < f < r and 
V(t) = Vi 7^ otherwise, which gives 



h{t) = IM) +lm{e- ieV ^- T ^ h J(t-T; V U V 2 )} 

-Im j e -<<*Wft e -ieV, J V X , V 2 ) } . (10) 

For t < 0, the corresponding current is obviously Idc(V±). 

The resulting 1(f) curves are depicted in Fig. [|. As 
is clearly seen, the combination of a bias and a mag- 
netic field produces nontrivial interference patterns with 
several underlying time scales. These include td, t± — 
2irh/\eV ± [iB9iH\, and either th = 2irh/ hb9iH or t v , 
depending on whether one is dealing with a pulsed bias 
potential or a pulsed magnetic field. In particular, for a 
pulsed bias potential and moderately large dc magnetic 
fields, there are current oscillations in the fall-off current, 
with a characteristic frequency of lu = (iB9iH/h. Such 
effects should be observable in ultra small quantum dots. 



h(t) = I dc {V 2 ) - Im {e~ leV ^ h J(t; V lt V a )} 



(9) 
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